Dissipative dynamics of a harmonically confined Bose-Einstein condensate 
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We study the dissipation of the centre of mass oscillation of a harmonically confined condensate in 
the presence of a disorder potential. An extension of the Harmonic Potential Theorem allows one to 
formulate the dynamics from the point of view of an oscillating disorder potential. This formulation 
leads to a rigorous result for the damping rate in the limit of weak disorder. 
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Trapped Bose gases provide an ideal setting for the 
study of nonequilibrium phenomena in a many-body sys- 
tem. Some examples include condensate formation fol- 
lowing a thermal quench FA, Q , collective excitations as a 
function of temperature [3J and the relaxation of highly 
nonequilibrium vortex states 0, [B[ . In these, and many 
other situations, the underlying superfluidity plays an es- 
sential role in determining the dynamical behaviour. 

In some recent experiments [6|48(j the dissipative dy- 
namics of a Bose condensate in the presence of a disorder 
potential was studied. This perturbing potential is the 
vehicle by which the collective centre of mass motion of 
the condensate is dissipated by means of internal exci- 
tations. The situation is analogous to the motion of an 
impurity through a superfluid where it is found that ex- 
citations can be produced above a critical velocity 
Here it is clear that the relative velocity of the impu- 
rity and superfluid is the relevant variable; the motion 
of a heavy impurity through a stationary superfluid or 
the flow of a superfluid past a stationary obstacle are 
physically equivalent. 

In this paper we demonstrate that a similar symmetry 
pertains to a Bose gas trapped in a harmonic potential. 
Harmonic confinement leads to an equivalence between 
the motion of the condensate through a disorder poten- 
tial that is at rest relative to the confining potential, and 
the harmonic motion of the disorder potential itself rel- 
ative to the trapping potential and the condensate. We 
exploit this equivalence to formulate a rigorous theory 
of the centre of mass motion in the presence of a disor- 
der potential and obtain an estimate of the damping in 
the limit of weak disorder. Our results are in qualitative 
agreement with those obtained earlier from an a naly sis 
of the one-dimensional Gross-Pitaevskii equation [lfjj. 

The Hamiltonian of the system studied experimentally 



H 



N 

E 

i=l 

Ho 



2m 



+ Vtrap(fi) 



Vw 



N 

E 

i=l 



Vdis(rj) 



(1) 



where Vtrap(r) = ^m(cj^p 2 +w^z 2 ) is the trapping poten- 
tial and Vdis(r) is a disorder potential whose properties 



we specify later. The interactions between the atoms is 
contained in Vint- To begin, we rephrase the Harmonic 
Potential Theorem (HPT) [11| in a form which will be of 
particular utility in the subsequent development. Start- 
ing with the many-body state |$), we define the state 
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|*) = T(x,p)|*), (2) 



where R = -4 Yli=i ^ s the centre of mass coordinate 
and P = J2iLi Pi is the total momentum of the par- 
ticles. These variables satisfy the commutation relation 
[Rp, P v ] — ihS^. The unitary operator T shifts the state 
in position space by x so that |^(ri, ...,i"Ar)| 2 = |$(ri — 
x, ...,rjv — x )| 2 - At the same time, the state is shifted 
in momentum space by p/N so that |^(pi, Pjv)| 2 = 
|$(pi — p/N, pN-p/N)\ 2 . This implies that the total 
momentum is boosted by p. 

We now consider the time evolution of the state \^>) 
according to the unperturbed Hamiltonian Hq. Using 
the Hcisenberg equations of motion for the operators R 
and P with respect to Hq, we find that 

|*(<)) = e - lMot/h \y) = f (x(t), p{t))e- lSot/h \^), (3) 



with 
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where M = mN is the total mass of the system. If 
we now take |$) to be an eigenstate |$ Q ) of Hq, we see 
that |*(n,...,rjv,t)| a = |$ a (n -x(t),...,rjv - x(i))| 2 . 
In other words, the probability density rigidly follows 
the motion of the centre of mass of the system. This 
is essentially the content of the HPT. More generally, 
the system can be described by the density matrix p — 
J2 a Pa\^a)(^a\- If the states |$ Q ) are all eigenstates of 
Hq, the total density of the system for the density matrix 
T(x, p) / 5T 1 '(x, p) oscillates rigidly according to n(r,t) = 
no(r — x(i)) where no(r) = Tr(pn(r)). This applies to 
the special case of a thermal equilibrium distribution. 

We next consider the dynamics of the centre of mass as 
governed by the full Hamiltonian including the disorder 
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potential. The Heisenberg equations of motion lead to 
the equation 
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where Z(t) = (V{t)\R z \^(t)) and F(t) = (* (t)\F z \^> '(*)) 
with F z = —J2iLidV<iis(*i)/dzi- Eq. © is an exact 
statement of the centre of mass dynamics, but requires 
knowledge of the dynamical state \^(t)). To determine 
this state we go to the interaction picture and define 
|*r(f)) = exp(iH t/H)\^/(t)) which satisfies 

!*,(*)) = |*(o)> -if dt> v di8 At')\yi(t')) (7) 




with V disJ (t) = exp(iH t/h)V d i s exp(-iH t/h). 

In the experiment we consider |7[, the centre of mass 
motion of the condensate is initiated by a sudden shift of 
the confining harmonic potential in the ^-direction. To 
describe this situation, we define the Hamiltonian of the 
system for t < to be 
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while for t > 0, the system evolves according to the 
Hamiltonian (QJ. The trap potential in ((SJ is illustrated 
by the dashed curve in Fig. QJa); we assume that the 
state of the system at t = is |*(0)) = |* ), the 
ground state of H' . This Hamiltonian can be expressed 
as H 1 = T(x, p)HT' ! (x, p), where x = zqZ, p = and H 
is 
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The external potentials of this Hamiltonian are illus- 
trated in Fig. QJb). The state |* ) = ft(x,p)|* ) is 
the ground state of H. 

Using the assumed initial state in J7|, we find that the 
state \^i(t)) — T^(x, p)|^j(t)) satisfies the equation 
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We thus see that is the state that evolves from 

| "Jo) a s a result of an oscillating disorder potential. 



FIG. 1: (a) The condensate, originally in equilibrium with the 
unshifted trap (dashed) and the disorder potential, begins to 
oscillate about the centre of the shifted trap (solid), (b) The 
condensate, originally in equilibrium with the trap and disor- 
der potential, is driven by an oscillating disorder potential. 



These results imply that the force appearing in ([5]) can 
be expressed as 

F(t) = <tt(t)|£|*(t)> = (h(t)\F,,j(t)\h(t)), (12) 



where 



F- , ( / i = -e iSat / H Y dV ^ i+x ^ e -iS t/H_ , , , 



Eq. (|T2|) is a key result and shows that there is an in- 
timate connection between the two very distinct physi- 
cal situations depicted in Fig. [1] In the first, one starts 
with an excited state corresponding to a displaced con- 
densate. This state then evolves according to J7]) in the 
presence of a static disorder potential. Even though the 
condensate follows a damped trajectory that eventually 
ends with the cloud being in equilibrium with the static 
disorder, the total energy of the system is conserved dur- 
ing this evolution. In the alternative situation described 
by (fl"0|) . the condensate starts in its ground state and is 
driven by a dynamic disorder potential moving according 
to the unperturbed centre of mass motion. In this case, 
the dynamic perturbation continually excites the conden- 
sate and the total energy increases as a function of time. 
However, the atomic cloud eventually reaches a steady 
state in which it moves with the disorder potential with 
no further increase in energy. That the cloud experiences 
the same force due to the disorder in these two situations 
is by no means obvious and is a consequence of the fact 
that the system is harmonically confined. 

More quantitatively, the solution of ((6]) for the assumed 
initial conditions is 
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where, because of (fl2|) . the second term on the right hand 
side is in fact the displacement for the state |^j(t)}. From 
(fT4]) we see that the change in the centre of mass position 
over one period T — 2tt/uj z is 



AZi = Z(Ti) — Z{Ti-\) 
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dt sinuj z tF(t), 
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where T; = IT. This is valid for both the static and 
dynamic disorder potential scenarios. 

To analyze the effects of F(t) on the dynamics of the 
centre of mass motion we will assume that the damping 
it gives rise to is weak. In keeping with this assumption, 
we evaluate F(t) perturbatively. To second order in the 
disorder potential we have 

F(t) = ($ |^,/WI*o> 



dt' e*' ($ \[F z>I (t),V diStI (*L,t')}\$ ) 



-~ / dt'^ \[F zJ (t),V d[s j(K(t')Xmo), (16) 



where |$o) is the ground state of Ho. The second term 
on the right hand side of (|16[) involving the positive in- 
finitesimal r] accounts for the lowest order effect of the 
disorder on the ground state l^o), while the third term 
arises from the dynamic perturbation of \^i(t)) in (fTUfl . 
We now write 

V disJ (x(t),t) = J drV dis (r + x(t))ni(r,t), (17) 

where hi{v,t) is the density operator in the interaction 
picture. The disorder potential is represented as 
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where the Fourier amplitudes U(k) are stochastic vari- 
ables having the following disorder averages: 



[/(k) = 0, U(k)U* (k') = (2tt) J <S(k - k')R(k). (19) 

Inserting (|17j) into (|16j) and performing the disorder av- 
erage, we find that F(t) = Fi(t) + F 2 {t) with 

Ftf) = ij dt'e^'J ^R(k)k z e^^x(K k; t ~ t>), 

FW) = if dt'J ^R(k)k z e*-W-^ X (K k; t - t% 

(20) 



where x(k, k;t—t') is the Fourier transform of the density 
response function 



X (r, r'; t - t') = -6{t - t 1 ) ($ 1 [hi(r, t), n/(r', t')]\$ ). 



The disorder averaged force in (|20|) is the main result of 
this paper and will be used to estimate the damping of 
the centre of mass motion in the linear response regime. 

In the experiments 0], the speckle pattern is one- 
dimensional so that R(k) = {2ir) 2 8{k x )8(k y )R(k z ) with 
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where a is the correlation length of the gaussian disorder 
and V^ 2 is is the square of the standard deviation of the dis- 
order potential. In this situation, the response function 
of interest is 

X (k z ,k z ;r) = [d 3 r [ d\'e lk ^- z "> X (r,v'-T). (23) 



A formal expression for this quantity can be given in 
terms of the exact Bogoliubov excitations of the system. 
However, here we make use of a local density approxi- 
mation (LDA) whereby each element along the length of 
the condensate is treated as part of a uniform cylindrical 
condensate having a density per unit length of v(z). In 
this approximation, the response function is taken to be 
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x(fc,fc;T)~ / dzx cy \{k,T;v{z)), 



where Xcyi(^> r ; v i z )) is the density response function of a 
uniform cylindrical condensate. Here and in the following 
we drop the z subscript on k z for convenience. In the 
Bogoliubov approximation, 
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where ipj(k) = 2-7T J Q dp p5rij{p,k) is the cross-sectional 
average of the mode density fluctuation Srij(p, k); the in- 
dex j distinguishes the various radial modes of the cylin- 
drical condensate. 

We calculate the density fluctuation by treating 
the condensate in the Thomas-Fermi (TF) approxima- 
tion 12J. This is a good approximation in the ex- 



perimental context since the number of atoms in the 
cloud is of order 10 6 . In the hydrodynamic limit, the 
normalization of the density fluctuation is then given 
by 2-7T J Q R± dp pSrij(p,k) — hiOj(k)/2g, where R±(z) — 
\\J R 2 Z — z 2 is the transverse TF radius at the position z 
along the axis. Here, A = u z /uj±, R z = ^2p/muj 2 , p is 
the chemical potential and g — Aitah 2 /m. 

It can be shown that the contribution of F\ (t) is negli- 
gible in comparison to F2(t). We thus focus on the latter 
in the following. Substituting ([24)) together with (|25|) 
into ([2"U| wc obtain 
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FIG. 2: The damping rate vs. vo/cq. The parameters used 
in this calculation are iT]: u z /2tt = 3.6 Hz; oj ± /2tt = 180 
Hz; the s-wave scattering length a = 200 ao; fi/h = 1 kHz; 
a = 10.6 /im. The solid curve is for the cylindrical LDA while 
the dashed curve is for the bulk LDA. 



For the initial conditions being considered, z(t) 
zq cos uj z t. In this case, we have 
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where w„ = ntu z and J n (x) is the Bessel function of the 
first kind of integral order n. This Fourier expansion 
is substituted into (|26p and the resulting expression for 



i*2 (t) is used in ([15} to evaluate AZi explicitly. Remark- 
ably, we find that AZi is virtually independent of I; the 
I = 1 result differs from the I — > oo limit by a few per- 
cent. There is essentially no transient on the time scale 
of T and implies that Z{Ti) ~ z (l - blT). Defining the 
damping of the oscillation as bT = —AZ^/zq, we thus 
find 



We now compare our results to the measurements pre- 
sented in Fig. 3 of Ref. 0]. The initial displacement 
of the harmonic potential of ^700 fim corresponds to 
Vq/cq — 2.9. For the weakest disorder strength given 
of V^Jli 2 = 0.0064 we find (b/u) z ) tb ~ 0.03, whereas 
(6/w z ) cxp ~ 0.04 [7J. This should be taken as reasonably 
good agreement given that there are no adjustable pa- 
rameters in the calculation. In this regard, we emphasize 
that the damping rate cannot be adequately character- 
ized using a white-noise spectrum 14j . 

We have also analyzed the data of Ref. Q which is ob- 
tained using a different protocol to excite the centre of 
mass oscillation. Here, the disorder is switched on sud- 
denly only after the oscillation of the condensate has been 
initiated. We find that the linear response damping rate 
is still given by ([28| in this case. Using the experimental 
parameters corresponding to Fig. 7 of Ref. 8], we obtain 
(6/w,)th - 0.06 for v a /c = 2.9 and V^J^i 2 = 0.0064, 
whereas (b/u) z ) exp ~ 0.002 [§]. We have no explanation 
for this discrepancy. Perhaps the analysis of a different 
situation such as a gaussian perturbation [8] may shed 
light on the limitations of the linear response calculation. 

In summary, we have shown that the dissipative dy- 
namics of the centre of mass motion can be formulated 
in terms of a conventional response function approach 
even though the initial state of the system is far out of 
equilibrium. With the appropriate response functions, 
this approach can also be used to study the dissipation 
at finite temperatures and in fermionic systems. 
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where Vq — uj z zq. 

The results of our calculation of the damping rate 
based on (|28|) are shown in Fig. [2j We observe that 
the damping rate exhibits a resonant peak at vo/cq ~ 1, 
where cq = y/ ///2m is the sound speed in a cylindrical 
condensate with chemical potential \i [IH • Also shown in 
the figure is the result obtained using the bulk LDA [l3[ 
whereby each element of the condensate is treated as a 
homogeneous gas. The cylindrical LDA is an improved 
approximation since it explicitly accounts for the effect 
of the transverse confinement on the excitations of the 
system and becomes exact in the limit of a small aspect 
ratio A (io z — > with /i held fixed). 
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